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ON STABILITY AND STABILIZATION OF ELASTIC 
SYSTEMS BY TIME- VARIANT FEEDBACK 



Abstract. We study a class of elastic systems described by a (hy- 
perbolic) partial differential equation. Our working example is the 
equation of a vibrating string subject to linear disturbance. The 
main goal is to establish conditions for stabilization and asymptotic 
stabilization by applying a fast oscillating control to the string. In 
the first situation studied we assume that system is subject to a 
damping force; next we consider the system without damping. We 
extend the tools of high-order averaging and of chronological calculus 
for studying stability of this distributed parameter system. 



1. Introduction 

Stability and stabilization by imposing fast oscillations onto a system have been 
extensively studied in literature. Basic example is stabilization of equilibrium of 
inverted pendulum by means of fast harmonic oscillation of its suspension point, 
PQ. Methods of high-order averaging developed in [2] are applicable to study 
stability and asymptotic stability for a wider class of (non-harmonically) time- 
variant systems. As an illustration, a stabilization condition for pendulum under 
non-harmonic (fast) oscillation of its suspension point has been established. An 
extension of this work onto the case of planar and spherical double inverted 
pendulum was done in [3]. 

In the present work we apply the latter approach to a class of distributed 
parameter systems: we study a possibility to stabilize an elastic system (a string) 
destabilized by a linear disturbance. A condition for asymptotic stability under 
application of a state feedback control to the damped system is presented. We 
also obtain a stability condition for the undamped vibrating string under an 
output feedback control. 

Our main results are Theorcm l5.1l in Scction[S]and Thcorcm lB.ll in Section 
El Our problem setting in Sectionals close to the one of 0] but the class of 
stabilizers and the results differ; we provide some comments on it in Section|7| 
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2. Notation 

Let SI C K be a closed interval. As usual L2(fl) is the space of all measurable 
(in the sense of Lesbegue) functions on Q having finite norm 

||u|| = u 2 dx 

The inner product in £2^) is denoted by < •, • >. The Hilbert space consisting 
of the elements £2^) which possess generalized derivatives up to order m also 
in £2(0), is denoted by Wf^O). The norm in this space is denoted by || • \\ m - 
Bellow, W^i^i) stands for a subspace of the space W^i^l) which coincides with 
the closure of the set of infinitely differentiable functions with compact supports 
in Q. 

Let {(j) n (x)} be a orthonormal basis of L 2 {Q). We denote by Gn be the 
subspace of £2^) spanned by {4>n(x)}n=i and by G N the orthogonal comple- 
ment of Gn- For an element v £ ^2(^2) we denote by vn its projection onto Gn 
and by v N the projection onto Gn- 




3. Problem setting 

We consider forced elastic system described by a (hyperbolic) partial differential 
equation. In particular, we consider the equation for an elastic string 

(1) u TT = a 2 u xx , x S O 

where denotes the partial derivative of u with respect to a variable £ and the 
domain = [0, 2ir}. We fix boundary conditions 

(2) u(x,r)=0, xedn. 

We introduce into JQ) a disturbance j 2 u and the equation becomes 

(3) u TT = a 2 u xx + 7 2 u 

whose zero solution is unstable if 7 2 > a 2 /A. 

If in addition, we assume that there is a damping, then Q turns into 

(4) u TT = a 2 u xx — au T + 7 2 it, x € f2 

with damping coefficient a > 0. 

We would like to achieve stabilization (asymptotic stabilization) of the 
null solution of © by applying time- variant high-gain linear feedback 

control term 

(5) h(u,kr) — 5 k 2 g(kr)u. 
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The resulting "controlled" equations is 

(6) u TT = a 2 u xx + 7 2 u + h(u, k r). 
and 

(7) u TT = a 2 u xx — au T + j 2 u + h(u,kr). 

The function g(r) is assumed to be bounded, continuous and 1-periodic, 5 > is 
a small parameter and k > 1 can be chosen arbitrarily. Besides, we will assume 
for g: 

Assumption 1. J* g(£) d£ = 0; 

Assumption 2. For G(t) = / * <7(£)^£ there holds 

/ G(£)d£ = 0. 
Jo 

Observe that, under these assumptions, the time-average of the vibra- 
tional term vanishes, in a way that the averaged equations © and (7J, coincide 
with the equations © and respectively. Therefore these averages are un- 
stable. 

For given initial data 

(8) u(x,0) = tp(x), u T (x,0) = ip(x), 

where ip € (fi)nW 2 1 (^) an d ^ G W2 (^) the corresponding (classical) solution 
for ©,© and ©,© exists and is unique in W 2 2 (Qt), <2t = x (0,T), 0. By 
a (classical) solution we mean an element of Wf (f2) which satisfies the equation 
for all £ and for almost all cc G fi. Next we precise the meaning of stability and 
asymptotic stability of the null solution. 

Definition 1 (Stability). The null solution u(x,t) = of ©,© is 
stable if for every 77 > there exists < £ < rj such that for any solution u(x, t) 
with \\u(x,0)\\l + \\u t (x, 0)|| 2 < £ implies 

||u(:M)|| 2 + |M:M)|| 2 <r] 

for alH > . ■ 

Definition 2 (Asymptotic stability). The null solution u(x,t) = 
of 10 j© is exponentially asymptotically stable if it is stable and in addition 
there exist <r, C > such that 

H-,t)|| 2 + K(.,i)|| 2 < Ce^WHxMll + lk(z,o)|| 2 ), 

for all t > 0. ■ 
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For convenience wc perform the change of the time variable t = kr. The 
transformed equation © is transformed into 

u tt = k~ 2 a 2 u xx + k~ 2 7 2 u + k~ 2 h(u, t) 
u(0,t) = u(2ir,t) = 0, 

while for the damped equation we obtain 

Uu — k a 2 u xx — k^ 1 aut + k~ 2 "f 2 u + fc~ 2 h(u, t), 
^ u(0,t)=u(2n,t) = 0. 

The initial data (jHJl is converted into 

(11) u(x,0) = <p(x), u t (x,0) = kr 1 ip(x). 

4. Chronological calculus: variational formula, logarithm and stabil- 
ity 

Chronological calculus has been developed by Agrachev and Gamkrelidze in |H] 
for studying nonlinear time-variant systems. The systems, to be considered in 
this work, are linear and we present some reformulations for the (time- variant) 
linear case. 

Consider the time- variant linear system of ordinary differential equations 

(12) z(t) = A(t) z(t), z{0) = z o 

where z(t) £ R™ and t ^ A(t) e R nx " is a matrix- valued function depending 
continuously on t. Considering its solutions z(t; Zq) we can introduce a flow of 
linear maps P : zq i— » z(t, Zo), or a fundamental matrix solution. Obviously, P* 
is the unique solution to the matrix ordinary differential equation with initial 
condition 

P l = A{t)P\ P° = Id. 

Following IS], we call the flow P 4 left chronological exponential of A(t) 
and denote it 



P* = exp* f A(t) dr. 
Jo 



We can also define right chronological exponential exp J Q A{t) dr as solution of 
P* = P t A{t), P° = Id. 

Using this notation, the solution of i|12|) is represented as 



z(t) = exp / A(t) dr z = P* z . 
Jo 



Stability of elastic systems 



■5 



The flow exp J* A(j) dr admits a Volterra series expansion of the form 

exp I A(r)dr = Id+ / A(r 1 )dT 1 + 
Jo Jo 

+ V / / . . . / A(n) . . . A(r n ) dr n ... d n . 
1 Jo Jo Jo 



4.1. Variational formula 

Assume now that, in (|12fl . we take A(t) = B(t) + C(t), considering B(t) as a 
reference matrix of coefficients and C (t) as its perturbation. Then the following 
chronological calculus variational formula holds, jS], 

(13) c5cp [ [B(t) + C(t)] dr = exp / B(r)dTo 

Jo Jo 

oe5q>^ adB(6)d6^J C(r)dT. 

Here the operator "ad" corresponds to the matrix commutator a.d A(t) B(t) = 
-[A(t), B(t)} = B(t) A(t) - A(t) B(t) and Q* = exp / * adB(6») d6» is the solution 
of the operator differential equation 

4<9* = Q t oadB(t), Q° = Id. 
dt 

4.2. Formal expansion for In exp L A T dr 

def 

Let A(t) — A t be a matrix-valued function with time-variant entries. We are 
interested in a formal expansion for the logarithm 

Ao, t(A T ) = In exp / A T d T . 

Jo 

It has been shown in [S] that Aq^(A t ) admits a representation 

00 

(14) A 0< t( A r) = E ^ 

m— 1 

where 

(15) A(™) = fdn fdT 2 ... f m 1 dT mPm (A Tl ,.. . A Tm ), 



and, for each m > 2, p m (A Tl , . . . ,A Tm ) is a homogeneous polynomial of first 
degree in each A Ti . Moreover, it is a commutator polynomial in A Tl , . . . , A Tm , 
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i.e., p m can be expressed as a linear combination of A Tl , . . . , A Tm and of their 
iterated commutators. 

Expressions for the first A^ m ^ can be founded e.g. in 0; in particular 
A^ 1 ) coincides with the averaging of A T : 

(16) A« = f A T dr. 

Jo 

The series 1)140 is known to be absolutely converg ent if / || Ar|| dr < 0.44, 0. 

4.3. The logarithm and stability of time-variant systems 

Consider linear fast-oscillating system 

(17) z(t)=A(kt)z(t) 

where x £ K™, A(t), t > 0, is matrix-valued function, which is continuous and 
1-periodic with respect to r; k > is a large parameter. 

A condition for stability of time-periodic system 117|) uses the monodromy 
matrix P 1 . If all the eigenvalues of P 1 are located in the interior of the unit 
circle, then the system is asymptotically stable; system is unstable if at least 
one eigenvalue lies outside unit the circle. In general it is difficult to compute 
spectrum of P . 

Due to the asymptotic expansion l|14|l it is more convenient to deal with 
the logarithm InP 1 . One can arrive to a conclusion on the system's stability 
by analyzing the truncations of the logarithm series if one gets, somehow, an 
estimate for the rest term of such truncation. 

In terms of logarithm, stability conditions can be formulated as follows: 
system (|17|) is asymptotically stable if all the eigenvalues of In P 1 are located in 
the open left complex half-plane and is unstable if at least one eigenvalue lies in 
the open right complex half-plane. System (fTTjl is stable if all the eigenvalues of 
InP 1 have non positive real part and purely imaginary eigenvalues are distinct. 

5. Vibrating string with damping 

We consider the problem of asymptotically stabilizing the null solution of l|10|) 
by means of a state feedback control 

h(u,t) = Sk 2 g(t)u. 

In other words we study asymptotic stability of the null solution of the equation 

u t t = k~ 2 a 2 u xx - k^ 1 au t + [fc~ 2 7 2 + 5 g(t)] u, 
( ' u(Q,t)=u(2ir,t) = Q. 
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5.1. Infinite dimensional system of ODE 

Following the Fourier method, we look for a solution of (|18|l of the form 

u(x,t) = T{t)X(x). 

We will use dot and prime, respectively, to denote the derivative with respect 
to t and to x. Differentiating u and replacing it and its derivatives in the partial 
differential equation we obtain, for X ^ 0, T ^ 0, 

f + k^af + b{t)T _ X" _ 

(19) kF r o?T ~^~ X ' 
where 

(20) b(t) = -[k- 2 1 2 + Sg(t)}. 

From (|19|) it is clear that A must be constant, that is X" — XX, X £ R and 

(21) f + k- 1 af+[b(t)~Xk- 2 a 2 }T = 0. 
The problem 

(22) X" = XX, X(0) = X{2 tt) = 0, 
is a "spectral problem for the elliptic operator" 

(23) ^ = £ 

whit boundary conditions derived from the original ones in ijTHJ. We conclude 
that in J22J) 

O 

A„ = -Mn = — J") ^(i) = sin(/x„x), 71 = 1,2,... 

are the eigenvalues and eigenfunctions. Instead of 4> n (x) we take 

(24) cj) n {x) = tt-^Mx); 

now {4>n} forms an orthonormal basis in £2(0). 

We denote by T n the solution of (|21|l with A replaced by X„ = —fJ^- The 

series 

00 

(25) u(x,t) = Y / Tn(t)4>n(x) 

n=l 

satisfies, formally, l|18l) . For given initial data we obtain 

00 00 

71 — 1 71 — 1 
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Since {<f> n } is an orthonormal basis in L2(£l), there holds 

(26) T n (0) = and T„(0) = fc" 1 (V, K)- 

Therefore, if T n (Q) and T„(0) are given by (|26[1 . then series (|25[1 formally satisfies 
JTSJ with given initial data. 

Substituting u, given by (|25|) into the partial differential equation l|18|) 
and equalizing the terms which contain sin(/i„ x) we obtain an infinite system 
of second order linear ordinary differential equations 

(27) f n (t) + k- 1 a f n (t) + [b(t) + ^4 k- 2 a 2 } T n (t) = 0, 

for n — 1,2,.... Formally speaking, i|18|) is equivalent to infinite-dimensional 
system of equations (I27|l . 

In the next section we illustrate the essence of our method by studying 
asymptotic stability of an equation (|27|l . 

5.2. Stabilization of (J23) 

If we take time- averaging of the coefficients of (|27fl . then, due to Assumption ^ 
the average of b(t) equals — k~ 2 7 2 and we arrive to the equation 

(28) f n (t) + k- 1 af n (t) + k- 2 (ii 2 ia 2 - ^T^t) =0. 
If destabilizing linear disturbance is sufficiently large, e.g., 

2 ^ 2 2 

7 > « Mn 

then the averaged equation (|28|) is unstable. We prove though that, even in this 
case, one can achieve stability by choosing sufficiently large "frequency" k of 
(non-harmonic) oscillation J5|. 

To derive an asymptotic stability condition for a single equation in i|27f) we 
first transform it into a two-dimensional first-order linear system by introducing 
a new variable S n — kT n (c.f. [3]). We get for = (T n , S n ) 

(29) z n (t) = A n (t) z n (t) 
where 



An(t) — , _! , 2 ,:?. ,,2 ;.- 1 . 



o fc- 1 

We represent A n (t) as B„(i) + C„, where 



B n (t) = 6kg(i)B, B 

and 





1 



— ™ I 2 2 2 

1 7 -a -a 
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According to the chronological calculus variational formula l|13(l , the mon- 
odromy matrix of the system (|29|l can be represented as 

P„ x = el? / A n (0 d£ = ex? / [B n (0 + C n ] 
Jo Jo 

= exp y B n (r) dr o exp ^ ^exp J ad £?„(#) d6 C n J dr. 

From the Assumption ^ it follows that 

/ B n {r) dr = O2 and, hence exp / B n (r) dr = ld.2, 

Jo Jo 

where O2 and M2 are, respectively, the two-dimensional null and identity ma- 
trices. Besides 

5xp aAB n {6)d6 \ C n = e - 5fcG « adB C„ 

where f Q g(d) dd = G(t) (c.f. Assumption EJ. 
We represent the monodromy matrix as 

P* =cxp / D n (t)dt, 
Jo 

where 

D n (t) =e - 5fcG (*) ads C„ 

= [Id 2 -<5 fc G(t) ad B + G 2 (i) ad 2 £ - ffl^ G(t) ad 3 B + ...] C n . 

As far as ad J SC„ = 2 , j > 3, the series for e - 5kG ^ adB C n ends at 
the third term and therefore 



(30) D n (t) = 



5G{t) k- 1 
fc- x (7 2 - a 2 /4) - G(t) -k^a-SGit) 



where G(t) = a5G{t) + k5 2 G 2 (t). 
Let now 



A = In Pi = In ^cxp J D n (t) dt 



be the logarithm of the monodromy matrix. 

Consider the expansion (jl4|l for the logarithm, with A^ m ^ defined by (|15|l . 
We will show that, under some conditions, the stability properties of are 
determined by A^. 
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As we said in Section f4.2l coincides with the averaging of D n (t). 
According to Assumption [21 

(31) A(1) = ( k-^^-JuD-SHT -k-^ct ) 
where 

(32) r= f G 2 {t)dt>{). 

Jo 

The eigenvalues of A^ can be either real or complex conjugated; they 
are equal to 



(33) e - 9±V T^ 

where = — fc _1 a and A = 8 2 V — k~ 2 ("f 2 — a 2 /i^) are the trace and the 
determinant of A^ 1 ), respectively. Obviously, 8 < 0. A^ possess real negative 
eigenvalues if 



(34) k-^-^^KS 2 [ 1 G 2 (t) d t<k~ 2a2+ ^ 2 - a2n2 . 

4 Jo 4 

For large n the latter inequality is never satisfied and then for 

(35) S 2 J^G 2 {t)dt>k- 2 a2+ ^ 2 - a2n2 

we obtain a pair of conjugate complex numbers with negative real part. 

What for the rest term A — A' 1 ', then the following estimate is available. 

Proposition 1 Q3 ). For n 2 < 4 fc + ^ there exists a positive constant 
c and a 2 x 2 matrix R such that 

(36) K-k {1) =S 2 R, \\R\\ M <c 
where neither c nor R depend on n. ■ 

Here || • ||m stands for some matrix norm. The upper bound on n is 
irrelevant; we can assume it to be fulfilled for k > fco as stated in the next 
proposition. The estimate (|36|1 can be derived from convergence results for the 
chronological exponential presented in [H]; in [3] the detailed exposition for the 
particular case of stability of inverse double pendulum can be found. 

From (|34|l . (|35l) and H36fl . we arrive to the condition for the asymptotic 
stability of the two-dimensional system (|29() . 
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Proposition 2. For each e > there exist So > 0, k > 1 such that 
the null solution z n (t) = of (|29|> is exponentially asymptotically stable if 
< 5 < Sq 7 k > k and 

(37) S 2 f 1 G 2 (t) dt > fc -2 4 7 2 - a 2 n 2 + ^ 



and unstable if < S < (5q, k > ko and 



(38) S 2 f G 2 {t)dt < k 

Jo 



1 /l„,2 „2^2 



_ 2 — a n 



4 



Observe that if the condition (|37() is fulfilled for n — N then it is fulfilled 
for all n> N. 



5.3. Asymptotic stability of Galerkin's approximation 

The Galerkin approximation for l|18|) is a system of (decoupled) equations l|27() 
for n= 1,...,AT. 

This system is equivalent to a 2 TV block diagonal system 

(39) z N {t) = H N {t)z N {t) 
where zjj(t) — (zi(t), . . . ,ZN(t)) and 

H N (t) = diag{D 1 (t),...,D N (t)}. 

Due to the decoupling, we may conclude that the system l|39|l is asymptotically 
stable whenever each of the systems 

z n (t) = D n (t)z n (t), n = l,...,N, 

is asymptotically stable; it is unstable if at least one of them is unstable. Insta- 
bility occurs if l|38|l holds for some n = n\ < N (then it holds for all n < ni). 
According to Proposition [2] system l|39() is exponentially asymptoticaly stable if 

(40) S 2 f G 2 {t)dt>k- 2Al2 ~ a2 1 - 



4 

Observe that the latter condition does not depend on N . 

Let Gn be the subspace of Z/2(^) spanned by {4>n(x)}n=i defined in 
and let be the orthogonal complement of Qn in L2{&). 

We denote the projection of the (unique) solution of l(T%|) on Qn by ujy 
and the projection of u onto Qj^ by u . 

Then 

TV N 

u N {x,t) = ^(u, <j) n )<j) n (x) = T n (t)4> n (x) 

is the solution of the Nth Galerkin's approximation of l|18l) and the following 
proposition holds. 
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Proposition 3. For each e > there exist S > 0, k > if < 6 < 
Sq , k > ko and 

(41) S 2 J^G 2 {t)dt>k- 2 ^ 2 ~ a2 +e 

then the null solution z n (t) = of the Galerkin's approximation of l|18|) is 
exponentially asymptotically stable I 

Therefore, if ijHJ) holds then 

\\u N (-,t)\\l + \\(u N ) t (-,t)\\ 2 ^0 as t^+oo. 

To establish asymptotic stability of the problem 0,101 it suffices to prove, that 
for sufficiently large N, the projection of u onto Qjj tends to zero as t — ► +oo 

\\u N (-,t)\\ 2 + ||uf(-,t)|| 2 ^0 as t^+co. 
5.4. Convergence of \\u N (■ , t)\\l + \\u^(-,t)\\ 2 

As far as u is an element of W^iQr), where Qt = f2 x (0,T), u N (-,t) is in 
W|(fi) for all i G (0, T). We prove now 

Proposition 4. For the projection u N there exists c > such that 

+ ||<M)|| 2 < e- c *(||^(-,0)||? + \\v» (,0)|| 2 ), 

for all t > 0. ■ 

Proof. Introduce a function y defined, for any u £ W 2 (Qt), by 

(42) V- = — ^ — 1|«|| 2 H ^ — || W;c || 2 + -|| Ui || 2 + — <U: u t > . 

We recall that < •, • > and || • || stand for the inner product and for the norm 
in Lzip), respectively (c.f. Section Let Vn denote the restriction of the 
function V onto Qh. 

We first note that there is a positive constant C\ such that 

(43) V N >C x {\u N f y + \\u¥f). 
Indeed, V can be represented as 

y = ^_o_||^||2 + i|| Ut ||2 + 1 Wk^au + utf . 

As long as the inequality ||u^|| 2 > -/V^u^H 2 holds, then 

v N > ^Vim 2 + ^n^ii 2 + l\\u?r + \ iifc- 1 *^ + uf f 
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and 63 follows with d = d{k) = min [^-, \j ■ 

Let Vjv denote the time derivative of Vn{u n (-, t)). After some simple 
computation we conclude 

-V N > C 2 (N, k,t)\\u N \\ 2 + ^ a 2 \\u?\\ 2 + ^IKI 

where 

C 2 (N,k,t) = ^a 2 N 2 - 

As long as g(t) and b(t) are periodic and bounded, there exists (3 > such that 
IK*) | ^ for a11 1 ^ °> ( ci - W- Then > for all t > 

C 2 (N,k,t) > C 2 (N,k) = ^-^a 2 N 2 - 

We assume N to be sufficiently large to guarantee C 2 (N, k) > (for all S < So, 
k > ko). Then, for those N, C 2 can be chosen independently on N 

-Vn > C 2 (N,k)\\u N \\ 2 + a 2 \\u N x \\ 2 + ^Ikfll > 0. 

We can also show that 

^<^-V|m 2 + ^a 2 ||^|| 2 + || U f||. 

Therefore, there exists a positive constant C(k) such that Vn < — C(fc)Vjv, 
hence 

(44) V^(*)<^(0)e- c W*. 
From (|4*H|l and we conclude 

\\u N {-,t)\\l + \\u*{;t)f as t^+oo. ■ 

From the Propositions|3]and0]we conclude the asymptotic stability of the 
null solution of infinite dimensional system 139|) and the consequent asymptotic 
stability of the null solution of the problem Q,©- 

5.5. Main result for damped string 

Sumarizing the reasoning of the previous sections we state 

Theorem 5.1. Consider the equations Q,© f° r damped elastic string 
subject to linear disturbance and controlled by a time-variant (high-gain fast- 
oscillating non- harmonic) linear feedback control Q. 

For each e > there exist So > 0, ko > such that if < S < So, k > ko 
and the inequality (|41() holds, then the null solution u(x, t) = of problem 
lj7)l.P]l is exponentially asymptotically stable. ■ 



i^K*) + 7^& 2 (*) 

2 k L a 
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6. Undamped vibrating string 

In this section, we consider the problem of stabilization of the null solution of 
undamped equation We try to stabilize it by a forcing term, which results 
in equation 

Uu = k~ 2 a 2 u xx + k~ 2 7 2 u + fc~ 2 h(u, t) 

(45) 

u(Q,t) = u(2n,t) = 0. 
Our choice is a time- variant "output" feedback, , 

(46) h(u,t) = 5k 2 g{kT)U N u, 

c.f. (O, where II at stands for the orthogonal projection of the "state" u onto 
the subspace of L2^l) generated by the first N harmonics. 

For the control h defined by (|46[) , the equation l)45|l is formally equivalent 
to an infinite dimensional system of second order linear ordinary differential 
equations 

(47) f n (t) + [k- 2 (rf l a 2 - 1 2 )-5g{t)}T n (t)=0, n<N 
f n {t) + k- 2 ( M 2 a 2 - 7 2 ) T„(t) =0, n > N. 

This system is obtained by introducing the (formal) expansion (|25(l of u. 

The analysis of the stability of un — TInu goes along with what was done 
in Sections 15.21 and 15.31 for the damped case. 

As before, we introduce a new variable S n — kT n and reduce l|47|) to a 
two-dimensional first-order linear system 

(48) z n {t) = A n (t) z n (t) 
where 

A(t)-( ° 

An[t) ~ \ k-itf-nla^+kSgit) 

Using chronological calculus variational formula (|13f) . we compute the mon- 
odromy matrix for (|48|l 



Qli = cxp / 
Jo 



D n {t)dt, 



where 



D(t)-( SG ® ) 

Un[t) -\k- 1 { 1 2 -a 2 tf l )-k5 2 G 2 {t) -8G(t) ) 



and G(t) — f Q g(£) d£. Next we consider the logarithm A = lnQ* given by the 
asymptotic expansion l|14H and compute A^ 

A (D -( o fc- 1 
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where V = f Q G 2 (t) dt > 0. It is clear that the main difference between damped 
and undamped systems is that now is a traceless matrix while (|31f) had 
negative trace. Hence, to achieve stability we have to ensure that its eigenvalues 
are two imaginary (conjugate) numbers. Stability condition for (each of) the 
two-dimensional system ()48|l is analogous to the one of the Proposition [21 

Proposition 5. For each e > there exist So > 0, k > 1 such that 
the null solution z n (t) = of (|48|) is exponentially asymptotically stable if 
< S < So, k > k and 

and unstable if < 5 < S , k > fc and 

S 2 [ 1 G*(t)dt<k-^ 2 - a2n2 - e . ■ 
Jo 4 

Recall that Galerkin's approximation for the equation (|45|l is a system 
of N decoupled second order equations |@BJ, where n = 1, • • • N. As in Section 
15.31 this system is stable whenever all the equations are stable, that is, when 
the Proposition [5] holds for all n < N. 

To be able to conclude stability of the null solution of (|45|l , 146|) we have 
to analyze stability of the projection of (|45|l onto Qjj 



Let us consider a function 

v = -k-^\\u\\ 2 +k-^\\u x r+ 1 -\\u t r. 

This function is not necessarily positive but its restriction Vm onto Qj^ is positive 
provided that N is sufficiently large, say, TV 2 > 2j 2 /a 2 . Indeed, since the 
inequality || zt^ || 2 > A^ 2 1| || 2 holds, then in this case 

v N >^ [« 2 ^v 2 - 7 2 ] \\u N f + fc-^KII 2 + i||«ril a > o. 

Let Vn denote the time derivative of Vn(u n (•, t)). Then 

V N = -fc- 2 7 2 < /,«f > +k- 2 a 2 < u»,u» t > + < u»,u% > 

= -k- 2 7 2 < u N ,u? > -k- 2 a 2 < u^ x ,uf > +k- 2 a 2 <u?,u» x > + 
+ fc- 2 7 2 <v?,u N > 
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The equality < , u^ t > = — < u xx , uf > , which holds under the boundary 
conditions is concluded by integration by parts with respect to x. 

Therefore Vn(u n (-, t)) = const wherefrom we obtain an upper estimate 
for u N and conclude stability of the null solution of of the projection of l|45|) onto 
Gjj. Completing this by the result of the Proposition we conclude stability of 



Theorem 6.1. Consider the equations for an undamped elastic 

string which is subject to linear disturbance and is controlled by a time- variant 
(high-gain fast-oscillating non-harmonic) linear output feedback control (|46|l . 
Then for each e > there exist Sq > 0, k > such that if < S < So, k > k 
and the inequality 



holds, then the null solution of problem ©,101 is stable. ■ 



7. Final remarks 

From (J2HJ, it is clear that the degree of asymptotic stability increases with a. 
An interesting but not unexpected phenomena arises when 2 — 4A equals zero: 
multiple negative eigenvalue bifurcates to a conjugate complex pair. 

If l|38|l holds for some n — m < N, then finite dimensional approximation 
(|39|l is unstable implying that the null solution of the problem 101, © is unstable. 

The output feedback control constructed in the Section is suitable for 
the asymptotic stabilization of the damped string (system 0,10)). 

The approach to the problem presented in this work is more general than 
the one of 0|. The main difference is that we do not select a particular type of 
vibrational control, in contrast to the harmonic vibration and rather restrictive 
assumption for frequencies, which appeared in 0]. Besides here we were able to 
deal with stabilization of an undamped string by means of output feedback. 
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